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Abstract 

We establish the equivalence of a recently introduced discrete model 
of the hydrophobic interaction, as well as its extension to continu- 
ous state variables, with the Ising model in a magnetic field with 
temperature-dependent strength. In order to capture the effect of 
symmetries of the solvent particles we introduce a generalized multi- 
state model. We solve this model - which is not of the Ising type - 
exactly in one dimension. Our findings suggest that a small increase 
in symmetry decreases the amplitude of the solvent-mediated part of 
the potential of mean force between solute particles and enhances the 
solubility in a very simple fashion. High symmetry decreases also the 
range of the attractive potential. This weakening of the hydrophobic 
effect observed in the model is in agreement with the notion that the 
effect is entropic in origin. 
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1 Introduction 



The hydrophobic interaction is a property of aqueous solutions which has a 
long and multi-faceted history of study [|l], ||. We focus here on one par- 
ticular issue, viz. the solvent-mediated interaction between solute particles. 
The balance between an energetically favorable accommodation of a solute 
particle between solvent particles and the entropically unfavorable necessity 
to place these solvent molecules in some particular orientation leads to an 
increase in free energy which is believed to be the basic mechanism of the 
hydrophobic effect. Since this increase of free energy becomes weaker with 
smaller separation of solute particles an effective solvent-mediated hydropho- 
bic attraction arises which is superimposed on and may even dominate the 
direct interaction between solute particles. 

This interplay between favorable energy and unfavorable entropy has been 
illustrated in an idealized lattice model where each solvent particle is allowed 
to be in one of q discrete states (orientations) 0] . There is a nearest neighbor 
interaction between solvent particles with energy w if both particles are in 
one particular state and u (with u> w) otherwise. A solute particle may be 
accommodated in the interstitial site between solvent particles, provided both 
solvent particles are in that special state, which is energetically favorable by 
the amount u—w. Calculation of the solvent-mediated part of the potential of 
mean force between solute particles and of the solubility (done exactly in one 
dimension 0] and numerically in two and three dimensions Q]) has revealed 
in all space dimensions an inverse relationship between the magnitude of the 
attractive force and its range. Furthermore, the solubility of the hydrophobe 
decreases with increasing temperature. In one dimension the decay of the 
potential of mean force with distance is monotonic, while in two and three 
dimensions there is an oscillatory modulation. 

Here we first review and extend this basic model to continuous state vari- 
ables, thus being somewhat more realistic with regard to the interpretation of 
the molecular states as orientations (Sec. II) and show that both the original 
discrete model and the continuous model are equivalent to the ferromagnetic 
Ising model with a magnetic field with a temperature-dependent strength 
(Sec. III). Since this mapping retains the local structure of the interactions, 
all the known results about the Ising model can be translated directly into 
properties of the hydrophobic model. Finally, we model symmetries of the 
solvent particles in a simple fashion by allowing for the accommodation of 
solute particles between solvent molecules that are in any n rather than just 
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one energetically favorable pair of special states (Sec. IV). The resulting 
(n + l)-state model is not equivalent to an Ising model, but may still be 
solved exactly in one dimension. These exact results are discussed in the 
concluding Sec. V. 

2 Model of the hydrophobic attraction with 
discrete and continuous state variables 

In the model of Refs. 0, |^ each solvent molecule is placed on a d-dimensional 
lattice and can be found in q different states. These states may be thought 
of representing orientations with respect to some arbitrarily chosen reference 
axis. There is a nearest neighbor interaction which favors pairs of molecules 
which are both in one particular state. We shall refer to this state as the 
special state, or state 1. Thus the interaction energy between molecules in 
states 2, j has the form 



Without loss of generality we shall set the global reference energy m = and 
therefore consider the case w < 0. The energy of interaction of an interstitial 
solute molecule with its neighbors is taken to be f . This quantity, however, 
does not enter the expression for the solvent-mediated part of the potential 
of mean force; only the solubility depends on f ||^. 

We calculate the potential of mean force W{r) on the basis of the potential- 
distribution theorem [^]. We define Pu as the probability of finding a given 
neighboring pair of molecules in the special state 1 and P{r) as the prob- 
ability of finding two such pairs at a distance r, measured in lattice units. 
Following one then finds 




w, when i = j = 1 
u, otherwise 



(1) 



with 



u > w. 



(2) 



W{r) = -kTln [P(r)/p2j 



(3) 



at temperature T. 

As solubility we define the dimensionless ratio 



S — Psoln/ Pf 



gas 



(4) 
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where psoin is the number density of the solute in the solution and Pgas is 
its number density in an ideal gaseous phase in osmotic equilibrium with the 
solution. For small solubility (such that the saturated solution is very dilute) 
this quantity is the Ostwald absorption coefficient. Within the lattice model 
one finds [§] 

S = Pne-"/'^. (5) 

The equilibrium properties of the model are encoded in the partition 
function 

^L = Ee-^''^''^ (6) 
c 

where L is the number of lattice sites, j3 = l/{kT), and the sum is over all 
configurations C the system may be found in. The nearest neighbor sum 

E{C) =wY: (7) 

<i,j> 

is the total energy of a configuration C of state variables Sj = 1, . . . , g. 

For the actual computation of the quantities Pn and P(r) one uses the 
standard transfer matrix formulation of statistical mechanics models; see e.g. 
0. In one dimension the transfer matrix is the q x q matrix 



(8) 



which has 

(9) 

at the (1,1) position and 1 everywhere else. The partition function becomes 

Zl = trace V^. (10) 

In order to calculate thermal averages we also introduce the projectors 
S{k) = I A;)( A; I on states k. These are the diagonal matrices with 1 at the 
diagonal element k and elsewhere. This yields the formal expressions 

Pn = ^trace[S(l)VS(l)V^-i] = A(i|v^-i|l) (n) 
P{r) = ^trace [S(l)VS(l)V'-i)S(l)VS(l)V^-"-i] 

= |^(l|V'-^|l)(l|V^-"-^|l) (12) 







1 ... 1\ 


V = 


1 


1 ... 1 




[ 1 


1 ... 1 ; 


b -- 


= e~ 


-w/kT ^ -j^ 
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In the thermodynamic hmit L — > oo the leading contribution to the par- 
tition function comes from the largest eigenvalue A of the transfer matrix. 
Anticipating the finite gap in the spectrum of V we may therefore write 

Pn = (13) 
P{r) = Y,Qu{r) (14) 

where pi is the probability of finding a molecule in state 1 and the correlation 
function 

Qn(r) = lim ^trace [S(l)V'-^S(l)V^-''+^] (15) 

is the probability of finding molecules at both sites (i,i + r — 1) in state 1, 
irrespective of the configuration of other molecules. 

In higher dimensions the transfer matrix is constructed analogously. On 
a hypercubic lattice in d-dimensions the transfer matrix can be written as a 
product of two non-commuting matrices V = V^^^V^^^. The diagonal matrix 

yd) = e-'^^ (16) 

gives the Boltzmann weight of a configuration in a (d — l)-dimensional hy- 
perplane with 

Y^^wY. 'S,(1)S,(1). (17) 

<ij> 

The projectors Sj act as unit matrices on all sites except i. The prime on 
the sum marks the restriction of the summation to the hyperplane of states 
on which the transfer matrix acts. The non-diagonal part 

v(^)=nvi (18) 

i 

is the tensor product of the one-dimensional transfer matrix over all sites 
of the hyperplane. Thermal averages are calculated by inserting projectors 
Sj(A;) at the appropriate positions in V and in the trace over V^. 

In the previous discussion the states that the molecule may be found in 
were chosen to belong to a discrete set. Since the model is not designed to 
be realistic for any particular solvent such an artificial representation of the 
orientation of a molecule is not really a cause for concern. Nevertheless, for a 
proper interpretation of these states as orientations it is useful to extend the 
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model to continuous state variables by the following limit procedure. One 
imagines there to be n distinct states 1, . . . , n which have mutual interaction 
energy w. Thus there are ri^ pairs of states with energy w and — pairs 
which have interaction energy m = 0. Hence we shall refer to this model as 
the n^-model. These n states may be interpreted as labeling infinitesimal 
solid angle segments on a sphere surrounding the solvent molecule (see Fig. |I| 
for a two-dimensional representation). Not all of these segments have to be 
in close proximity on the sphere; they may belong to distinct blocks which 
define those areas on the surface of the molecule between which a solute 
particle may be accommodated. Taking the limit n, g — oo with the ratio 
q = q/n fixed, one obtains a continuous version of the model where the range 
of states [0, 47r/ g] represents the total amount of solid angle on the surfaces of 
each of the two neighboring solvent molecules between which solute molecules 
may be found. Parts of that interval may then refer to distinct areas on the 
surface of the solvent particle. 

For fixed n and q the one-dimensional transfer matrix has a block form 
similar to (^), but where instead of the matrix element h one has an n x n 
matrix with all elements equal to h. All other matrix elements are equal to 1. 
As will be shown in the next section, the n^-model and the original 1^-model 
are equivalent, thus proving that the earlier restriction to discrete states did 
not constitute an additional oversimplification. In the interpretation of the 
1^-model as a model with a continuous degree of freedom the special state 1 
simply represents a fraction 47r/g of solid angle (oriented in some direction) 
rather than a discrete single direction. Notice, however, that the analogue of 
the quantity Pn is the double sum 

^11 - E P^, (19) 

which turns into an integral in the limit n, g oo. The pair correlation 
function P(r) has to be redefined accordingly. With the proper redefinitions 
all results of the 1^-model may be interpreted in terms of the continuum 
limit of the n^-model, i.e., with the integer q being replaced by g = g/n, and 
1 < g < oo considered a real number. 
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3 Equivalence to the Ising model 



In the n^-model all special states are indistinguishable, as are the non-special 
states. Hence it is important to know only to which group of states the state 
of a molecule belongs, but not to which particular state within each group. 
Any joint probability involving some state k depends only on whether k 
belongs to the group of special states {k G {1, . . . ,n}) or not {k E {n + 
1, . . . ,q}). Hence the n^-model is equivalent to a 2-state model, i.e., to an 
Ising model. We may identify the special states with an Ising spin Si = +1 
whereas the other states may be collectively identified with Sj = —1. 

In order to work out the equivalence in detail we note that the partition 
function (^) of the original 1^-model may be cast in the form of a double sum 

^e-mc)^^ ^ (20) 

C Ci C'iCi) 

The first sum runs over all fixed configurations Ci of A^^i molecules in state 1; 
the second sum then includes all other configurations C with the configura- 
tion Ci specified by the first sum. Since the energy E{C) = E{Ci) depends 
only on the positions of molecules in state 1, i.e., on Ci, the second sum can 
be evaluated trivially to give 

E =(g-l)"-"^^ (21) 

C'(Ci) 

This originates in the degeneracy of states with a fixed number A^^i of molecules 
in state 1 in the configuration Ci, and we obtain 

Zl = Ee-^^'^^^n? - l)^""^^ = Ee-^^^""^^ (22) 

Ci Ci 

with a new energy function 

E = E + kTln{q-l)Y,{l- 5,^,i). (23) 



In those terms the partition function (p2D of the model for hydropho- 
bic attraction turns into the partition of a ferromagnetic Ising model in a 
magnetic field with energy 

E = — J ^ SiSj — h'^Si — const. (24) 

<i,j> i 



7 



The Ising interaction parameters are given by the relations 

J = -J (26) 

h = -^--A:Tln(g-l), (26) 

with c the coordination number of the lattice. This enters since part of 
the contribution of the magnetic field strength is contained in the nearest 
neighbor sum (J^). If some lattice sites have a coordination number different 
from c (e.g. at the boundary in case of non-periodic boundary conditions), 
additional local magnetic fields would contribute to the Ising energy. The 
temperature dependence of the entropic contribution kT\n{q — l)/2 to the 
field arises since the number of unfavorable states does not depend on tem- 
perature and hence the factor kT must cancel in the partition function. The 
constant term in ( p^ does not affect the thermodynamical properties of the 
system and may be dropped. Notice that g is a discrete variable and therefore 
within this formulation of the model h changes discontinuously as a function 
of q. 

A similar resummation may be performed for the partition function of the 
n^-model by first fixing a configuration of favorable states with Nf molecules 
(which yields a factor (g — n)^~'^f from the summation over the remaining 
degenerate configurations) and then counting the degeneracy of all such fa- 
vorable states. For each given A^^i = Nf this yields another factor n^^ . Thus 
we are left with an Ising partition function with a magnetic field 

C1I1 1 

h' = -—--kT\n{q/n-l). (27) 

The continuum limit of the state variable may now be taken and yields an 
Ising partition function identical to that derived from the original 1^-model 
but where now q is is replaced by q = q/n and treated as a continuous 
variable in the range 1 < g < oo. Hence the n^-models are all equivalent 
to each other and to the Ising model in a magnetic field given by (|27|), i.e., 
by the replacement q — > q/n. Without loss of generality one may therefore 
discuss only the 1^-model. 

On the level of the transfer matrix description the equivalence may be 
shown by a similarity transformation V = YVY^^. Because of the cyclic 
property of the trace the partition function is invariant under any such trans- 
formation. For the 1^-model we define 

u = e2-/(^-i) (28) 
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The matrix 



-1 

1 

1 

1 




1 





1 



-2 




1 

a;-2(<?-2) 



V 



-ig-2) ^-2(g-2) 



LU 



-(9-2)2 



(29) 



has as its inverse the transposed matrix, = . The transformation 
leaves the projector S(l) invariant and brings V into a block diagonal form 
with the 2x2 matrix 

W^(^f^) (30) 

in the upper left diagonal and everywhere else. Therefore the trace reduces 
to a trace over the two-dimensional subspace on which V acts in a nontriv- 
ial way, i.e., trace V'^ = trace V'^ = trace W'^. In W one recognizes the 
transfer matrix of the Ising model and hence the equivalence is recovered. In 
higher dimensions the transformation Y is applied to each lattice site, thus 
reducing the full model to a two-state Ising model with parameters as defined 
above. 

Varying the temperature T corresponds to moving along a curve in the 
field-temperature plane of the Ising model which starts at T = 0, /i = 
— cw/4 > (all molecules are in the favorable state) and crosses the line 
= at a point defined by the condition 



w 



-kT ln(g - 1) 



(31) 



or, equivalently, b^^'^ = q — 1 (Fig|^). In two and three dimensions the Ising 
model has a phase transition of second order at finite temperature along this 
zero-field line. Above the critical point the probability of finding a molecule 
in the favorable state is 1/2, below the critical point one has spontaneous 
symmetry breaking between a "favorable" and an "unfavorable" phase. On 
a d-dimensional square lattice (c = 2d) some values of q at the critical point 
are given by 



1 + e 



-clw/kTc 




1 (exact) 

2 (exact) 

3 (numerical) 
oo (mean field) 



(32) 
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At the critical point the range of the potential of mean force would diverge 
if the model were applicable at such low values of q. At large distances the 
decay would be algebraic and proportional to the decay of the magnetization 
correlation function. 

For larger values of q (which seem more realistic) the zero-field condition 
corresponds to subcritical temperatures. Hence, following the curve defined 
by varying the temperature of the hydrophobic model one expects a jump 
discontinuity in the probability of finding a favorable state from a value 
greater than 1/2 to a value less than 1/2 at h = 0. The strength of the 
discontinuity increases with increasing q. The regime which appears to be 
relevant for the hydrophobic effect 0, ^ corresponds to h < 0. For all 
coordination numbers (c = 2 being the smallest possible) this is consistent 
with the requirement that the change of free energy AF = w + kT In (g — 1) 
arising from the accommodation of a solute particle be positive. 



4 Symmetry effects 

Implicit in the construction of the 1^-model is the view that the surface struc- 
ture of the solvent molecule is made up of certain regions between which a 
solute molecule may be accommodated. No provision is made for more struc- 
tured solvent molecules which have chemically different surface regions. One 
could imagine n non-identical regions on each solvent molecule such that each 
identical pair of regions allows for the accommodation of a solute particle be- 
tween two solvent molecules, but the accommodation between distinct pairs 
would be energetically so unfavorable that it may be considered forbidden. In 
this section we extend the model to allow for n pairs of orientations between 
which a solute particle may be accommodated (Fig. |T|). We denote these 
special states by 1, . . . , n, with 1 < n < q and assign non-positive interaction 
energies 

{w, when 1 < i = j < n . , 

0, otherwise 

We shall refer to this model as the n-model. With regard to solubility a 
solvent molecule of type n described by these interactions would have a higher 
symmetry than a molecule described by the original 1-model with the same 
value of q, since there are more, but mutually exclusive pairs of favorable 
orientations. The limiting case n = q corresponds to the zero-field, g-state 
Potts model We stress that we define this model in the spirit of Refs. 
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1^, ^, i.e., not with the intention of describing a particular real solvent, but 
with the aim of exhibiting what may be found to emerge as generic behavior 
resulting from an increase in symmetry. Since some of the important features 
of the 1-model were already seen in the simple one-dimensional model, we 
restrict our analysis of the n-model also to the one-dimensional case. This 
can be solved exactly by diagonalizing the transfer matrix. A reader not 
interested in the details of the diagonalization may skip the following more 
technical subsection. 



4.1 Exact solution of the one-dimensional model 

The q X q transfer matrix V of the n-model has b for its first n diagonal 
elements and 1 everywhere else. By a transformation similar to (|29|) - with 
the 1 at the (l,l)-position replaced by an x n unit matrix and the quantity 
g — 1 in the normalization factor and in ui replaced by {q — n) - the transfer 
matrix becomes block diagonal with 



W 



1 y/q — n \ 
1 y/q-n 



1 1 



b i/g — n 
^q — n q — n j 



(34) 



in the upper left corner and elsewhere. Hence the model is equivalent to an 
(n + l)-state model, which is, however, not equivalent to some Ising or Potts 
model. Nevertheless the degree of degeneracy is high and therefore it is not 
difficult to guess the eigenvectors of W and hence of V by a second trans- 
formation, again analogous to Y, but with the first two columns replaced 
by matrix elements of the form Yn = y^/z+, Y^i = and Y12 = y^/z+, 

= for all 2 < < + 1, and with the quantity q — 1 replaced by n. 
Since the columns of the transformation matrix are eigenvectors of W the 
variables 



y± 



2Vq 



(q + l~2n-b± ^J{b-q-ly + 4:n{b-l)] . (35) 

/ 71/ \ / 



are determined by the two solutions of a quadratic equation. The constants 



-± 



(36) 
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are normalization factors for the eigenvectors. 

The full transformation which diagonalizes V may then be constructed 
and all quantities of interest may be calculated. One finds four distinct 
eigenvalues which take the values 6 — 1 (which is (n — l)-fold degenerate), 
{{q — n — l)-fold degenerate), and 



The largest eigenvalue is A+. Before proceeding further we give the asymp- 
totic expansion up to first order in 1/g which is necessary for the study of 
the hydrophobic regime of large values oi q ^ b. One finds: 



The exact components A^'^^ of the eigenvectors | A ) belonging to the eigen- 
values A are given as follows: 




(37) 



A_ 



A 



q+ {b - l)n/q 
(6-l)(l-n/g) 
g-2(6- l)(l-n/g) 
[q + 2(6 — l)n/q] n/ {q — n). 



(38) 
(39) 
(40) 
(41) 



A = A±: 



y±/{z±^/q - n) {n + l<k<q) 



l/z± (1 < A; < n) 



(42) 
(43) 



A = 6 — 1: {p — 1, . . . ,n 



A?) 
A?) 




(44) 
(45) 



{n+l<k<q) 



A = 0: (p — 1, . . . ,q — n 



1; a;2 = e^'^^/^^-")) 



XiJ^) = (l<A;<n) 

Af = a;f^'-"-'V\/Q^ {n+l<k<q) 



(46) 
(47) 
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The matrix elements needed for the calculation of average values pre- 
sented below are then given by 



{k\V''\ k') 



4±E^ + 4^ (1< A; < n; n + 1< A;' < g) . (48) 
+ -4f^ (n + l< k,k' < q) 



These expressions are obtained by inserting a complete set of eigenstates of 
the transfer matrix, 

(A;|V"| A;') =^(A;| A)(A|V"| A;'), (49) 

A 

and using the components A*-'^-' = ( A; | A ) and their complex conjugates A*-^-** = 
( A I A; ) given above. 

4.2 Potential of mean force and solubility 

To study the properties of the model we consider first 

P{k) = lim trace (S(A;)V^) /Zl = \im {k |V^| k)/ZL (50) 

L—*oo L—*oo 

which is the probability of finding a molecule in state k. Using the matrix 
elements calculated in the above one finds 

s _Jl/2+ l<A;<n (special states) , . 

^ \ y+/izl{q-n)) n + l<k<q 

At the temperature defined by 

- w = A;T In (g - 2n + 1) (52) 

the probability 

n 

P, = J2P(k)=n/zl (53) 

k=l 

of finding a molecule in any one of the special states equals 1/2. This is the 
analogue of the zero field condition (^) in the 1-model in one dimension and 
suggests that the condition q — 2n + 1 > b limits the parameter range of the 
n-model relevant for modeling the hydrophobic effect. 
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The quantity entering the solubihty is the nearest-neighbor correlation 
function 



Pii = hm V trace {S{k)VS{k)y^-^)/ZL = — P,. 



(54) 



fe=i 



For the computation of the potential of mean force we also need 

n n 

Pir) = lim VV trace [S(A;)VS(A;)V'-'SmVS(nV^-'-VZi 
= lim X]E trace [P{k)Y'-^P{l)Y^-'-^]/ Zi 



■iL-2 



AlZ 



r-l 



The potential of mean force therefore takes the form 

W{r) = -kT In [l + yle-("-^)/« 

with the amplitude 



and the localization length 



In 



(55) 

(56) 
(57) 

(58) 



We remark that the expression of the potential of mean force does not involve 
the eigenvalue and eigenvectors of the the {n — l)-fold degenerate eigenvalue 
6—1 characteristic for the n-model. The eigenvalues X± and the quantities 
z±, however, do depend on n. 

The amplitude A is a monotonically increasing function of temperature, 
with A {Q/n — l)n^/6^ as T ^ and A = q/n — 1 infinite tempera- 
ture. The localization length vanishes in both these limiting cases and has a 
maximum at 6 = g + 1, independently of n. For this temperature we find in 
terms of the ratio u = n/q 



A 



1 + 



(59) 
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lni± 

1 - 



1/V4z/ 



(60) 

(61) 



For large q the maximal localization length grows proportionally to the square 
root of q, but decreases with increasing n. As a function of n the amplitude 
of the potential of mean force also decreases. The solubility is enhanced by 
a factor of n. 

In analogy to the original 1-model one expects the regime most relevant 
to the hydrophobic effect to correspond to q ^ b. Using again u = n/q and 
expanding the relevant quantities X±, up to first order in 1/g, (^) - (^1]) 
yield 



A 



V 



V 



In 



hv 



(6-l)(l 

-vl{kT) 



The inverse localization length may be written as a sum of two terms, 

h , 1 



i 



-1 



In 



In 



(62) 
(63) 
(64) 

(65) 



the second of which is positive and describes the reduction of as a function 
of n. As in the previous special case the amplitude decreases with increasing 
n, while the solubility increases by a factor of n. 

At sufficiently large distance the potential of mean force decays exponen- 
tially 

W{r) ~ -kTA&-^'-^^l^ (66) 

For small n and large q the localization length ^ depends on n only very 
weakly. The amplitude is reduced by a factor of n while the solubility is 
increased by the same factor. Hence one finds for two different systems the 
relations 

W{r) _ v' 

V 



W'{r) 



- = — for nv+w 



for ^ = ^' 

T'{v + w) 



(67) 
(68) 



15 



1/^ - 1/^' = In {q/q') for Tw' = T'w. (69) 

One can analyze in a similar fashion the effect of symmetry within the 1^- 
model where higher symmetry corresponds to a larger solid angle of accom- 
modation and hence to a decrease in q. This corresponds to setting u = 1/q 
in Eqs. - 

5 Concluding remarks 

The equivalence of the hydrophobic interaction model to the Ising model in 
a magnetic field with temperature-dependent strength allows for the inter- 
pretation of the solubility, the potential of mean force and other quantities of 
interest in terms of the well-known properties of the Ising model. For realistic 
values for water in the hydrophobic model, the strength of the magnetic field 
in the corresponding Ising model is significant, and consequently the model 
is far away from the Ising critical point. We can however not rule out that 
some real system is appropriately described by an Ising model, in which case 
a solvation effect associated with the corresponding short range interaction 
would change dramatically at the critical point. Indeed, the interpretation 
of g as a fraction of the surface area between which solute particles may be 
accommodated suggests lower realistic values of q than previously assumed 
and hence the zero-field condition for criticality might be satisfied if such a 
system exists. A more tangible consequence of a lower value of q (but not as 
low as required for the second order phase transition) would be to bring the 
minimum of the solubility as a function of temperature to larger and more 
realistic values, thus indicating that the mechanism by which this minimum 
occurs in a real system is similar to that described by the model. 

The most striking features of the original hydrophobic interaction model, 
the inverse relationship between amplitude and range of the potential of 
mean force and the decrease of solubility with increasing temperature in 
reasonable temperature ranges, are present also in the more symmetric n- 
model. However, the amplitude of the potential decreases with increasing 
degree of symmetry. This is in agreement with the notion that the origin 
of the hydrophobic effect is entropic. With increasing degree of degeneracy 
(and hence symmetry) the entropy decrease associated with the requirement 
to orient molecules in a particular direction is less and hence also the balance 
in the free energy difference is reduced. For solvents differing in their degree 
of symmetry the n-model predicts an inverse relationship (|67D also between 
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the ratios of the potential of mean force (at sufficiently large distance and 
temperatures at which the locahzation length are equal) and the ratios of 
solubilities (at temperatures such that the exponential prefactors are equal). 
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Figure 1: Accommodation of solute particles (black disks) on the interstitial 
sites between solvent particles (circles) on a square lattice. Depending on 
the model the figure has different interpretations, (i) In the n^-model with 
q discrete orientations (states) the whole circle consists of 4g discrete points 
which consist of quadruples related to each other by a 90 degree rotation. 
Each quadruple defines one orientation and there is always one quadruple 
of points which is located on the lattice axis. The special states l,...,n 
correspond to quadruples of points on the dark sections of the circumference 
of the circles. A solute particle may be accommodated whenever special 
points of two neighboring molecules are located on the lattice axis. (In the 
1^-model of Ref. there would be only one such quadruple of points.) (ii) 
In the continuum limit the interpretation of the open circles in the figure 
is that if dark segments of the circumference of neighboring circles inter- 
sect with the lattice axis then a solute particle may be accommodated, (iii) 
In the n-model defined below one distinguishes between n different special 
states. Accommodation of solute particles is possible only between neigh- 
boring solvent molecules which are in the same special state corresponding 
to one particular of n special quadruples of points. In the continuum version 
of the ra-model this corresponds to distinguishing between different surface 
areas on the solvent molecule. E.g. for n = 2 one would say that a solute 
particle could be accommodated if on both neighboring molecules either the 
long or the short dark segment intersect with the lattice axis, but not if for 
one molecule the long segments intersects and for the other molecule the 
short segment. (This scenario would allow for accommodation of a solute 
particle in the continuum version of the n^-model.) 
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Figure 2: Curve in the h — T plane of the Ising model which corresponds 
to varying T in the hydrophobic model (dashed line). Larger q leads to a 
steeper (more negative) slope. The bold line at h = marks the region of 
spontaneous symmetry breaking, ending in the critical point at Tg. In the 
regime h < relevant for the hydrophobic effect the probability Ps of finding 
a solvent molecule in the special state is less than 1/2. 
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